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Abstract — This paper solves reliability constrained unit 
commitment problem (UCP) for a composite power system 
network using network equivalent technique. Here, the 
integration of generation and transmission system reliability 
forms a composite power system model. The inclusion of load 
forecast uncertainty for the solution of unit commitment 
problem will give more accurate assessment of spinning 
reserve (SR). Probabilistic techniques help in setting the 
reserve requirement which is defined by reliability indices 
such as loss of load probability (LOLP) and expected energy 
not supplied (EENS). Reliability network equivalent 
techniques represent each market player in a power market. 
The unit commitment problem is solved by a genetic algorithm 
(GA) resulting in near-optimal unit commitment solutions 
and the required spinning reserve capacity is effectively 
scheduled according to the desired reliability level. The 
proposed enhanced inherited crossover operation in GA will 
inherit more information from the parent chromosomes 
thereby it improves the convergence fact and quality of 
solution. The effectiveness of the proposed technique for UCP 
is validated on IEEE RTS 24 bus system and a South Indian 
(Tamilnadu) 86 bus system. 

Index Terms — Composite generation and transmission system, 
Genetic algorithm, Load forecast uncertainty, Network 
equivalent technique, Probabilistic techniques, Unit 
commitment problem. 

Nomenclature 

Bj An outage condition of transmission network. 

C sti Total start- up cost of unit i during at k th hour. 
D k Total system demand at k th hour. 

DR(i) Ramp down rate limit of unit i. 
EENS 8 k EENS due to generation outage 
EENS kk EENS at each BLP at k ,h hour. 
EENS't EENS due to transmission outage 
EENS c k Composite System EENS at k th hour. 

EENS c k (s) Composite System EENS for the load step s of k th 

hour. 

EENS caltk System reliability index EENS at k th hour. 

F (Bj) Frequency of occurrence of outage B 

Fj(Pj, k ) Fuel cost function i.e. ^(I^)=a+b*I^ +c*P iJc 2 

©2011 ACEEE 
DOI:01.IJCOM.02.03.510 



a, b and c are generator cost coefficients 
/,,* Status of unit T at k th hour k 

H Dispatch period in hours. 

L Load curtailed due to generator contingency i. 

LC Total number of contingencies leading to load 

curtailment. 



"bj 



Load curtailed at BLP 'b' due to transmission line 



contingency j. 

LOLP 8 k LOLP due to generation outage. 

LOLP bk LOLP at each BLP 'b' at k ,h hour. 

LOLP'k LOLP due to transmission outage. 

LOLP c k Composite System LOLP at k th hour 

LOLP c k (s) Composite System LOLP for the load step s 

LOLP caIk System LOLP at kth hour. 

Pi State probability for state i. 

Piji Generation power output of unit i at k th hour. 

Pi,min Minimum power output of unit i. 

Pi.max Maximum power output of unit i. 

Pbf-bt Power flow between bus bf and bt. 

Pbf -bt, max Maximum power flow limit between bus bf-bt. 

P c j Probability of load that exceed the maximum load 

that can be supplied without failure 

P (B ) Probability of existence of outage B 

PL(s) Probability of load step s. 

SR k System spinning reserve at k th hour. 

T on (i) Minimum On time for unit i. 

T off (i) Minimum Off time for unit i. 

UR(i) Ramp up rate limit of unit i. 

V Pim in , V p max Minimum/Minimum limit of voltage at bus p. 

V p Voltage at bus p. 

X on (i, k) Time duration for which unit i is On at time k. 

X off (i, k) Time duration for which unit i is Off at time k. 
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I. Introduction 

Unit commitment (UC) plays a major role in the daily 
operation planning of power systems. System operators need 
to perform many UC studies, in order to economically assess 
the spinning reserve capacity required to operate the system 
as securely as possible. The objective of the UC problem is 
the minimization of the total operating cost of the generating 
units during the scheduling horizon, subject to system and 
unit constraints. The solution of the above problem is a very 
complicated procedure, since it implies the simultaneous 
solution of two sub problems: the mixed-integer nonlinear 
programming problem of determining the on/off state of the 
generating units for every hour of the dispatch period and 
the quadratic programming problem of dispatching the 
forecasted load among them. The evaluation of the system 
spinning reserve is usually based on deterministic criteria. 
According to the most common deterministic criteria, the 
reserve should be at least equal to the capacity of the largest 
unit, or to a specific percentage of the hourly system load. 
The basic disadvantage of the deterministic approach is that 
it does not reflect the stochastic nature of the system 
components. On the contrary, the probabilistic methods [1]- 
[4] can provide a more realistic evaluation of the reserve 
requirements by incorporating various system uncertainties, 
such as the availability of the generating units, the outages 
of the transmission system, and the load forecast uncertainty. 
These methods combine deterministic criteria with 
probabilistic indexes, in order to find a UC solution that 
provides an acceptable level of reliability. 

Most of the existing probabilistic methods [2]-[4], are 
based on the priority list (PL) method in order to solve the 
UC problem. The PL method is very simple and fast, but it 
results in suboptimal solutions. Dynamic programming (DP) 
[1] and Lagrangian relaxation (LR) [5] have also been used 
for the solution of the UC problem. The main disadvantage 
of DP is that it suffers from the "curse of dimensionality," i.e., 
the computational requirements grow rapidly with the system 
size. The Lagrangian relaxation method provides a fast 
solution but it may suffer from numerical convergence and 
solution quality problems. Aside from the above methods, 
there is another class of numerical techniques applied to UC 
problem. Specifically, there are artificial Neural Network, 
simulated annealing (SA), and GA. These methods can 
accommodate more complicated constraints and are claimed 
to have better solution quality. SA is a powerful, general- 
purpose stochastic optimization technique, which can 
theoretically converge asymptotically to a global optimum 
solution with probability 1 . One main drawback, however, of 
SA is that it takes much CPU time to find the near-global 
minimum. GA is a global optimization method which works 
well and efficiently on objective functions which are complex 
in terms of the nonlinearities and constraints imposed. One 
of the disadvantages of GA is premature convergence, 
because when the selection is based on the quality of the 
individual, the genetic information of the best individuals 
tends to dominate the genetic characteristics of the 
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population [6]. To avoid disadvantages of this type in UCP, 
in this paper, the use of a GA is proposed, which is based on 
an enhanced inherited crossover operation. Thereby, the 
exchange of genetic information between two individuals with 
the modified crossover operator is more and improves the 
quality of solution. 

The network equivalent technique is used to by integrate 
unavailability of generating units/transmission lines and load 
forecast uncertainty for the solution UCP. The generation 
system in the market is represented by an equivalent 
multistate generation provider (EMGP) and with an equivalent 
assisting unit approach [7] where reliability of generation 
system is calculated. Using the network equivalent 
technique, the transmission system is modeled [8] and 
inclusion of transmission system reliability is employed by 
computing the transmission line outage at each load point 
using conditional probability approach [7-8]. Then composite 
system LOLP and EENS can be estimated by adding the 
individual LOLP and EENS caused by generating system/ 
transmission system, respectively. The above reliability 
parameters are taken into account in order to provide a more 
realistic evaluation of the system reserve requirement. 

II. UC PROBLEM FORMULATION 

In this section, UCP is formulated by incorporating 
reliability constraints. 

A. Objective function 

The objective of the UCP is the minimization of the total 
operating cost of the generating units during the scheduling 
horizon. For N generators, the operation cost (OC) is defined 
mathematically as given in equation (1). 

H N N 

min OC = £ £ F t (P iJt ) * 7 W + £ C „ , ( j ) 

k=\ i=l i=\ 

B. Constraints 

The minimization of the objective function is subject to a 
number of system and unit constraints, as follows. 
1. Power balance constraints 

The total generated power at each hour must be equal to the 
load of the corresponding hour. This constraint is given by 



2. Spinning reserve constraints: 

N 

^P t , mia *Ii, k >D k + SR k ,kell,H] 

;=i 

3. Capacity limits of generating unit 

p < p < p 

^/,min — r i,k — M,max 

4. Unit minimum ON/OFF durations 

[X on (i,k)-T on (i)]*[Ii, k -i-Ii, k ]>Q 



r"ff 



[X°" (i, k-\)- T°" (/)] * [L k _i - L k ] > 



(2) 

(3) 

(4) 

(5) 
(6) 
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5. Unit ramp constraints 

Pi,k-Pi,k-i^UR(i) 

Pi,k-l-Pi,k^DR(i) 

6. Networks constraints 

a. Voltage magnitude constraints 

V ■ <V D <V 

y p ' min — ' P — y y?,max 

b. Line flow constraints 



\Pbf-bt ^ Pbf-bt,mnx , bf - bt < 



(7) 
(8) 

(9) 

NB (10) 

7. Reliability constraints 

Spinning reserve requirements in equation (3) can be 
calculated using either deterministic criteria or probabilistic 
techniques. However, in the proposed technique probabilistic 
reserve assessment with spinning reserve requirements are 
assessed according to the desired level of reliability specified 
for the composite generation and transmission network. 
Therefore spinning reserve requirement should satisfy the 
reliability constraints given as follows: 

-41H] (11) 



LOLP k < LOLP spec k 



EENS k < EENS spec k e [l,H] 



III. Reliability modelling 



(12) 



A. Modeling of Generating systems 

For the purpose of reliability analysis, each generating unit 
is represented by the two-state model [7], according to which 
a unit is either available or unavailable for generation. 
In [7], the notion of unavailability is referred to as the steady 
state outage replacement rate or, more commonly, as the forced 
outage rate. Consider a power system with T Generating 
system (GS) and therefore T EMGPs. Let 'm' be the number 
of BLPs. Each of the EMGP is connected to different BLPs 
through the transmission system. The capacity outage 
probability table (COPT) of an EMGP can be realized by an 
available capacity probability table (ACPT) [8]. The 
parameters of an ACPT are the state available generation 
capacity AG., state probability p and frequency f All EMGP 
(j=/l, . • ., h-1, h+/l, . . ., Y) are added one by one as new unit 
to generate Total system available capacity probability table 
(TSACPT). 
Calculation of Reliability indices 

The system reliability indices for generation outages are 
calculated at hour 'k' using the equations given below. 

LOLP s k =Y/) ; , ke[l,H] 



EENS g k = 'YjPjLj , k e [l,H] 



(MWh) 



(13) 



(14) 



A reduction in the time requirements of the proposed method 
can be achieved by omitting the outage levels for which the 
cumulative probabilities are less than a predefined limit [8], 
e.g., 10 7 

B. Modeling of Transmission systems 

Electrical energy is delivered from an EMGP to its Bulk 
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load points (BLPs) through the transmission network. Failures 
of transmission network components can affect the capacity 
that can be transferred from an EMGP to its BLPs. In the 
composite power system, it is virtually impossible to provide 
the same level of service reliability to every customer. There 
are two basic types of transmission failures which affect 
reliability level of transmission system at each BLPs. The 
first type of failure is caused due to the isolation between 
BLPs and EMGP. In this case, only the isolated bulk load 
points are affected and energy delivered at that BLP is 0. The 
second type of failure is due to the line flow constraint 
violations such as congestion. In this case, ac power flow 
has to be analyzed to determine network violations for each 
area. If there are no network violations, the energy that can 
be delivered at the BLPs can be up to their peak loads. If 
network violations exist, the loads for all of the BLPs have to 
be shed proportionally (due to equality transmission right) 
to release the network violations (corrective action) and the 
energy delivered for an isolated BLP will be its peak load 
minus the load shed [8]. Network equivalent for each 
transmission network for the corresponding BLP and EMGPs 
are shown in Fig 1. Equivalent multi-state transmission 
provider (EMTP) between EMGPs and BLPs is represented 
as EMTPs. Hence the total transmission system reliability is 
the average of the individual transmission reliability (EMTP) 
at each of the BLPs. 
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Fig. 1. Transmission system and EMGP 

Calculation of Reliability indices 

The conditional probability approach [7] is applied to 
these equivalent transmission models and reliability indices 
at each of the BLPs are calculated using equations (15)-(16). 

LOLP bk = ^[P(B ; )*/J 7 ] , b^[l,m] 
j=no. of states 



(15) 



EENS b , k = ^JP(B ] )*P ci *^MnM {MWh) 

}=naoj states 

System reliability indices due to transmission unavailability 
are calculated using following equations 



LOLP' k = max(]T WLP hk ) ke [1,//] 
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EENS'k = ]T(££A« W ), k € [l,tf] (MWh) 

6=1 

In this approach, only the single line outage is taken into 
account and the only assumption that is made during the 
calculation is that adequate generation is available and 
therefore loss of load at any of the BLPs is due to transmission 
unavailability (EMTP) between the corresponding EMGPs 
and BLP. 

C. Modeling of Composite system 

Composite system reliability index is calculated by adding 
the generation and transmission system reliability indices 
using equation (19)-(20). 



lation. When the crossover and mutation processes are com- 
( 1 °) plete a new population is generated which replaces the origi- 
nal population. This must be repeated until one of the con- 
vergence criteria defined for the problem is met. Each of these 
cycles is known as a Generation. 

A. Enhanced Inherited Crossover operation 

Crossover is an extremely important operator for the GA. 
It is responsible for the structure recombination (information 
exchange between mating chromosomes) and the 
convergence seed of the GA and is usually applied with high 
probability (0.6-0.9). The chromosomes of the two parents 
selected are combined to form new chromosomes that inherit 
segments of information stored in parent chromosomes. 



LOLP c b =LOLP s b +LOLP' b 



(19) 
(20) 



EENS c b =EENS g b + EENS' b 
D. Load Forecast Uncertainty model 

The short-term forecast of the load demand is a problem 
with forecast uncertainties that plays a major role in the 
solution of the UCP. The load forecast uncertainty is 
represented by a normal distribution curve in which the mean 
value is the forecasted peak load. The normal distribution 
curve is divided into a discrete number of intervals and the 
load in the midpoint of each interval represents the probability 
for that interval. To simplify the calculation, a seven-step 
distribution (0,±1<t,±2<t,±3<t) is often used where, standard 
deviation, is 2-5% of the expected load [7]. Such a model will 
encompass more than 99% of load uncertainty. Based on the 
above analysis, the LOLP and EENS indices for each hour 
can be calculated by the following equations: 

7 

LOLP calk = £ WLP c k (s) * PL(s), ke[l,H] (21) 

i=l 

7 

EENS alik =Y J EENg : k( S )*PL(s),kB[l,H] (MWh) (22) 



IV. Overview of Genetic Algorithm 

GA is general purpose optimization algorithm based on 
the mechanics of natural selection and genetics. They oper- 
ate on string structures (chromosomes), typically a concat- 
enated list of binary digits representing a coding of the con- 
trol parameters (phenotype) of a given problem. Chromosomes 
themselves are composed of genes. The criterion which evalu- 
ates the quality of each chromosome is given by the Fitness 
corresponding to the evaluation of each individual for the 
objective function. Once the fitness of each of the individu- 
als in the population is known, it is subjected to a Selection 
process in which the best evaluated individuals have a greater 
probability of being chosen as Parents for the exchange of 
genetic information called Crossover. Then a percentage of 
the Offsprings (individuals generated in the crossover) are 
subjected to the Mutation process in which a random change 
is generated in the chromosome. This mutation process pro- 
vides greater diversity between the individuals in the popu 
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Fig. 2 Traditional crossover operation 
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Fig. 3 Enhanced inherited crossover operation 





t 


S-tturt 

1 




) 






in gfcj xn l/_jj -—r n.± Lea 
p ial el* L ej T3 im e3 ict 


r = 


■3— L 



3 



1 p ' 






it ■■=!= 


= £ 


= 


~J 


v 


En 


Te 


— 


=■ -=p 


hen 
-=x- n/t 


t«3 


4- 



Baie'dein th-= F" itr. 








P L-= k Tap fix- xt p -= la x- -zj aaa -zj x -=a xxa -=■ x far 






t= =:-; t r a:ia =t =. t i=> ia 






+ 




IvJ ■— ana -=j i- Lae ■— x bent "v n. liac ■*.■>. -it h -cea ncip eaxa eliaa t£ 


= h ™=° ™*=M 




Fig. 4 Flowchart for proposed enhanced inherited GA 
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Traditionally, in GA a linear crossover is performed on the 
two chromosome represented as a string, as shown in Fig. 2. 
Here a chromosome 1 is combined with chromosome 2 of 
randomly selected chromosomes and two new chromosomes 
are formed from the information stored in parent 
chromosomes. Similar new chromosomes are formed between 
3 and 4 and so on. Here chance of inheriting information from 
parent chromosomes is less, since there may be chance to 
inherit more information while combining the chromosomes 
2 and 3 instead of chromosomes 1 and 2. In order to avoid 
this disadvantage, an enhanced inherited crossover operation 
is carried to transfer more information from parent to 
chromosomes. 

Fig.3 clearly shows the enhanced inherited crossover 
operation. At the end of the enhanced crossover operation, 
the number of population in GA is doubled, thereby the 
proposed enhanced inherited crossover operation inheriting 
more operation from the parent chromosomes. Then a 
percentage of the offspring (individuals generated in the 
crossover) are subjected to the Mutation process in which a 
random change is generated in the chromosome. This 
mutation process is carried for 2p size population and 
provides greater diversity between the individuals in the 
population when compared to simple genetic algorithm. When 
the crossover and mutation processes are complete, for new 
2p size population, the fitness is evaluated for each individual 
for the objective function. Once the fitness of each of the 
individuals in the population is evaluated, then individuals 
are arranged in descending based on the fitness values. Then 
first p population is picked up for the next 
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Fig. 5 Flowchart of the proposed reliability constrained UCP 

generation and process is repeated until it termination 
criteria is met. Here p is a population size. The flow of 
algorithm is shown in Fig. 4 
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V. Solution methodology 

The flowchart of the proposed reliability constrained UCP 
using proposed GA is shown in Fig. 5. 

VI. Numerical results 

Simulations are carried out on an IEEE RTS 24 bus system 
and a South Indian 86 bus system. All simulations are 
performed using a PC with Pentium(R) 3.40 GHz, 1GB RAM 
system on the MATLAB 7 environment. Three different 
strategies have been carried out to schedule the system 
spinning reserve and are given as follows: 

Strategy 1: Calculation of required spinning reserve by 
deterministic approach. 

Strategy 2: Calculation of required spinning reserve only 
with generation unavailability. 

Strategy 3: Calculation of required spinning reserve 
by considering generation / transmission unavailability 
and load forecast uncertainty. 

A. IEEE RTS 24 bus system 

The proposed method has been applied to IEEE RTS 24 
bus system which is also called as Roy Billinton test system 
(RBTS) and is given in [9] . The generation cost coefficient is 
adapted from [9]. Generation/transmission line reliability data 
are taken from [ 1 0] and load data is adapted from [11]. 
Strategy 1: 

In this section, unit commitment problem is solved using 
proposed GA algorithm and results are compared with [11- 
13]. Here the system spinning reserve is deterministically set 
to the maximum capacity of the largest committed unit as in 
the ref [11-13]. For the sake of comparison with [11-13], 
network constraints and reliability constraints are not 
included. Here UC problem is solved using the proposed GA 
and lambda-iteration technique is used to solve the economic 
dispatch problem. 

Table 1 gives the comparison of the minimum total 
operating cost obtained by GA with respect to other 
techniques reported in the literature. It is clearly seen that 
the proposed method yields better results than other 
technique so far proposed in the literature. The minimum 
cost so far reported in literature is $ 721352.9 [13] which is $ 
25.68 higher than that obtained from proposed GA. Also the 
proposed GA yields the better solution with the introduction 
of enhanced inherited crossover operation in GA. The UC 
status for case 1 is given in Table 2. 

In this section, reliability constraints are included along 
with all constraints. Here reliability constrained UC problem 
is solved using the proposed GA and ED problem is solved 
using conventional method. 

Table 1 comparison of results-strategy 1 



Solution Tschniqua-s 


Minimum valus (S) 


ANN-DP[ 1 1 ] 


723326.5 


ILRH2] 


725996.9 


IPL-ALH[ 13] 


721352.9 


Simple GA 


72716S.41 


EnaancBd Inli*rit=d crossover GA 


721327.21 
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Table II. uc schedule- strategy 1 



Hour 


Unit Status 1.2 26 


Lead 


Fad 

cost S 


St. CD st 


1 


11111111111100000000000000 


1700 


19015 


106.93 


2 


11111111111 100000000000000 


1730 


19371 





3 


1 1 1 1 1 1 11 110 100000000000000 


1690 


:s"S9 





4 


1 1 1 1 1 1 11 110 100000 100000000 


1700 


1B966 





5 


11111111111 100000000000000 


1750 


19610 


74.329 


6 


11111111111110000000000000 


1B50 


211B7 


130.53 


7 


11111111111111000100001000 


2000 


23E05 


172.02 


B 


11111111111111111000000000 


2430 


32643 


104B.3 


9 


11111111111111111000000000 


2540 


34726 





10 


11111111111111111111000000 


2600 


36050 





11 


11111111111111111111101110 


2670 


3B251 


115.37 


12 


11111111111111111111000000 


2590 


35B16 





13 


11111111111111111111000000 


2590 


35B16 





14 


11111111111111111000000000 


2550 


34917 





15 


11111111111111111111110000 


2620 


365B2 





16 


11111111111111111111100110 


2650 


<SjS 


74.5 B7 


17 


11111111111111111000000000 


2550 


34917 





IS 


11111111111111111000000000 


2530 


34534 





19 


11111111111111111000000000 


2500 


33963 





20 


11111111111111111000000000 


2550 


34917 





21 


11111111111111111111000000 


2600 


16050 





22 


11111111111111111000000000 


24B0 


335B4 





23 


11111111111111100100110000 


2200 


27423 





24 


11111111111110000000000000 


1S40 


21066 







TotalcDstfS)=7213. 


27.21 







Strategy 2: 

For sake of comparison with Ref. [14], loss of load 
probability (LOLP) index is used to set the level of spinning 
reserve. Also LOLP index is expressed in % and the lead time 
of the system is fixed for four hours. Also the network security 
constraints are not considered. Table 3 gives the comparison 
of the total operating cost obtained by proposed GA with 
respect to other techniques. The LOLP index is kept as 1.5 % 
for comparison. 

It can be inferred from the results that the solution 
obtained from the proposed method is better than the SA 
and GA. The proposed enhanced inherited crossover GA 
finds a solution which is $ 654 lesser for the LOLP 1.5 % 

spec 

when compared with SA and is found to be promising. The 
UC status for each unit of the 26 unit system for LOLP = 

J spec 

1.5 % is given in Table 4. 
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Strategy 3: 

In this section, three different case studies have been 
carried out to understand the importance of inclusion 
transmission line outage and load forecast uncertainty. 
Case 1 : Calculation of required spinning reserve only with 
generation unavailability. 

Table III. comparison of results- strategy 2 



E olution. T Khniqu ? 


Total operating cost S 


Minimum 

value 
710696 


Mean 

value 

710971 


Maximum 

value 

711315 


SA[14] 


Simpb GA 


710652 


712014 


713374 


Modified GA 


710042 


7106^ 


7110BB 



Table IV. uc schedule- strategy 2 



Hour 


Unit status L2....26 


Load 


Fuel 
Cost(S) 


St Cost 


1 


11111111100 00 0000000000000 


1700 


1B429 





2 


11111111100000000000000000 


1730 


1B7B5I 





3 


11111111000000000000000000 


1690 


:s:o: 





4 


1 1 1 1 1 1 11010000000000000000 


1700 


1B429 


S 1.606 


5 


liiiiinoioooooooooooooooo 


1750 


19030 





6 


11111111111000000000000000 


:s:o 


2047 B 


172.16 


7 


11111111111 100000000000000 


2000 


22B93 


132.62 


B 


11111111111111110000000000 


2430 


32071 


967.37 


9 


11111111111111110000000000 


2540 


3429 B 





10 


11111111111111111000000000 


2600 


35960 


360.62 


11 


11111111111111111000000000 


2670 


3760 B 





12 


11111111111111111000000000 


2590 


35726 





13 


11111111111111111000000000 


2590 


35726 





14 


11111111111111101000000000 


2550 


34541 





15 


11111111111111101111110000 


2620 


36345 





16 


11111111111111101111111000 


26:0 


37231 


39.993 


17 


11111111111111101000000000 


2550 


34541 





IB 


11111111111111110000000000 


2530 


34063 


27B.69 


19 


11111111111111110000000000 


2500 


33410 





20 


11111111111111110000000000 


2550 


34533 





21 


11111111111111110111000000 


2600 


35B00 





22 


11111111111111110000000000 


24 B0 


33025 





23 


11111111111111000000000000 


2200 


26617 





24 


11111111110000000000000000 


.540 


20262 







Total cost = 710042 


.64 S 







Table V. comparision of different casestudies (load =1700 MW)-strategy 3 





Nd 


DeterministicApproacli 


Probabilistic Approach 






10 % Df Load as reserve 


EENSspec=3 MWh 


EEN55pec=5MWh 




Cssel 


Case 2 


Case 3 


Casel 


Case 2 


Case 3 


Op. Cost (S) 


IB221 


18429 


1B449 


1B553 


IB659 


1B343 


1B449 


1B449 


LOLPcal 


16.99 


0.1B47 


0.1S46S93 


0.1S35631 


0.124659 


0.4600rB 


0.1B46S93 


0.1S46B93 


EENScal MWh 


127.6 


2.993726 


2.9934153 


1.6193416 


0.3924653 


4.6154566 


2.9934153 


2.9934153 


SR, MW 




222 


222 


29 B 


374 


146 


222 


222 


EENSspk. MWh 






j 


3 


3 


: 


5 


^ 
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Table VI. generated power schedule (load =1700 MW) - strategy 3 



Unii no. 


No R?s?rv? 


13 Mofloadas 5R 


EENSsp^= 3 MWh 


EENSspM.= 5 MWh 


Case] 


Cas-2 


Cas-3 


Case] 


Cas-2 


Cas-3 


Pe 


Re 


Pe 


Re 


Pi 


Re 


Pi 


Re 


Pi 


Re 


Pi 


Re 


Pi 


Re 


PI 

P2 
P3 
P4 
P5 
P6 
P7 
PB 
PS 
P10 
Pll 
P12-P26 


400 
400 
350 
145 
140 
135 
130 












400 
400 
350 
137 
132 
127 
123 
15.2 


15.2 











17.5 
22.7 

27.7 
32.5 
&>* 



60. B 






400 
400 
350 
137 
132 
127 
124 
15.2 
15.2 













18.1 
23 2 
27.6 
314 
60. B 
60. B 








400 
400 
349 
133 
128 
124 
120 
15.2 
15.2 
15.2 










1.27 

21.7 
26.7 
31 1 

34.8 
60. B 
60. B 
60. B 





400 
400 
342 
131 
126 
122 
11B 
15.2 
15.2 
15.2 
152 






7.92 
23.9 
2B.9 
33 2 
36.9 
60. B 
60. B 
60. B 
60. B 




400 
400 
350 
141 
136 
131 
127 
15.2 















14.2 
194 
23.9 
27.7 
60. B 










400 
400 
350 
137 
132 
127 
124 
15.2 
15.2 













IB 1 
23 2 
27.6 
314 
60. B 
60. B 








400 
400 
350 
137 
132 
127 
124 
152 
15 2 














1B.1 

23.2 
27.6 
31.4 

60. B 
60. B 









Case 2: Calculation of required spinning reserve by 
considering both generation and transmission unavailability. 
Case 3: Calculation of required spinning reserve by 
considering generation / transmission unavailability and 
load forecast uncertainty. 

In this study, only the EENS criterion is used to set the 
level of spinning reserve. In power system operation such a 
criterion is more appropriate than LOLP since EENS accounts 
for both the probability of outages and corresponding 
average load lost. For a system load of 1700 MW, Table 5 
compares the results using deterministic approach, with and 
without reserve and with two different EENS limit values 

spec 

in the probabilistic approach. It is observed from Table 5 that 
the cost of operation is well correlated with calculated 

level of spinning reserve (SR) where, the spinning reserve 
(SR) is found to increase as the EENS criterion (EENS ) is 
restricted. That is the SR , in case 3 is 222 MW when EENS 

cal spec 

is 5 MWh. Also SR cal - is increased to 374 MW when the 
EENS is 3 MWh. " 

spec 

From Tables 5 and 6, it is observed that the cost and the 
number of committed units are increased when transmission 
reliability constraints and load forecast uncertainty are 
included. Table 6 clearly gives the generation and reserve 
dispatch for each unit. When EENS is 3 MWh, an increase 

r spec 

in cost of 104 $/h is observed (Table 5) when moving from 
Case 1 to Case 2 due to increase in committed units(Table 6). 
Increase in cost signifies the importance of considering 
transmission line reliability issues for calculating spinning 
reserve requirement in power system model. Also when 
EENS is 3 MWh, an increase in cost of 106 $/h is observed 

spec 

(Table 5) when moving from Case 2 to Case 3 due to increase 
in committed units (Table 5). Therefore it is vital to include 
load forecast uncertainty when reliability issues are 
considered. It is to be noted (Table 5) that the EENS cal is well 
within EENS , in the probabilistic approach for all the three 
case studies. When there is no spinning reserve constraint 
it is obvious that the system is more unreliable with a LOLP 
of 16.999 and EENS , of 127.6 MWh. In deterministic method, 

cal 

the spinning reserve is considered as 10 % of the load. When 
SRis 10 % of the load, the actual levels of LOLP is 0. 1846 and 
EENS cal is 2.99 MWh. It is interesting to compare the 
deterministic reserve case when EENS is 5 MWh set for 

spec 

the probabilistic reserve case where both of them yield 222 
MW spinning reserve. When EENS is having different 

©2011 ACEEE 
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values i.e 3 MWh and 5 MWh, the reliability level for 
probabilistic approach will be higher or lower when compared 
with deterministic approach. However, increase or decrease 
in the reliability indices values depends on the choice of 
EENS . The choice of EENS can be determined by 

spec spec J 

probabilistic approach from the trade off between cost and 
reliability level which makes probabilistic approach 
advantages over deterministic approach. 

The commitment schedule for Case 3 when EENS is 3 

spec 

MWh is given in Table 7. Fig. 6 illustrates the variation of the 
available spinning reserve capacity over the scheduled period 
of 24 hours for the EENS is 3 MWh and 5 MWh. Obviously 

spec J 

the 3 MWh requires more spinning reserve when compared 
to 5 MWh. The total operating cost for various values of 
EENS limits are shown in Table 8. Here the total operating 

spec L o 

cost of the system increases as EENS criterion is restricted. 

Table VII. generated power schedule (24 hour) - strategy 3-case 3 



Hour 


Unitnc 1,2.... 26 


Load 


Fuel 
CostS 


St. 
lost S 


1 


11111111110000000000000000 


1700 


1B534 





2 


11111111110000000000000000 


1730 


1SB93 





3 


11111111110000000000000000 


1690 


'.S4'.: 





4 


11111111110000000000000000 


1700 


1B534 





5 


11111111110000000000000000 


1750 


19133 





6 


1 1 1 1 1 1 11 HI 100000000000000 


1S:j 


20 817 


221.0B 


7 


11111111111110000000000000 


2000 


23214 


132.62 


B 


11111111111111110000000000 


2430 


32071 


s. : . : .:: 


9 


11111111111111110000000000 


2540 


3429 B 





10 


11111111111111111000000000 


2600 


35960 


360.62 


11 


11111111111101111000000000 


2670 


37 BBS 





12 


11111111111101111000000000 


2590 


35992 





13 


11111111111111111000000000 


2590 


35726 


97.543 


14 


11111111111111101000000000 


2550 


34541 





15 


11111111111111101111100000 


2620 


3631: 





16 


11111111111111101111111000 


2653 


37231 


39.993 


17 


11111111111111101000000000 


2550 


34541 





IB 


11111111111111110000000000 


2530 


34063 


:"S.53 


19 


11111111111111110000000000 


2500 


33410 





20 


11111111111111110000000000 


2550 


34533 





21 


11111111111111110110000000 


2600 


35770 





22 


11111111 111 111110000000000 


24 B0 


33025 





23 


11111111111110100000000000 


2200 


269 S3 





24 


11111111111100000000000000 


1S40 


2069: 







Total operating co st = 7 


12546. 6t 
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B. South Indian (Tamilnadu) 86 bus system 

The proposed model is validated on South Indian 
TamilNadu 86 bus system. The generation data and network 
data is taken from [15-16] and the reliability data of generation 
and transmission system are taken from [17-18]. Due to page 
limitation data's are not included in this paper. The single line 

veACEEE 



diagram of a South Indian 86 bus system is given in Appen- 
dix A. The generating system is divided into four indepen- 
dent generating companies EMGP ( 1 to 4) . Using the network 
equivalent technique, the generation and transmission net- 
work is modeled and the simulation results for the various 
values of EENS limits are given in Table 9. It is noted that 

spec ° 

the spinning reserve required at each hour is increased or 
decreased in proportion to the specified reliability level. The 
commitment schedule for Case 3 when EENS is 50MWh is 

spec 

given in Table 10. 

Table VIII. total operating cost 





EENE,- =C =3 MWh 


EEXS, -,==:" MWh 




Cas-3 


Cas-3 


Co<S dav) 


"12:44.4: 


709922.93 
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4c: 



EC 




1 



tEN5ipec=5MWh 
E£N5ipec=3MWh 



3 5 7 » 11 13 15 17 19 21 23 

Hour 



Fig: 6 Spinning reserve at each hour for different EENS 
Table IX. total operating cost 





EENS,- ,=50 MWh 


EEN5,- = =10MWh 




Cas-3 


Cas-3 


CostCLNRdav) 


14870493 


149460^4 



Table X. dispatch schedules for case 3 (eensspec=50 MWH) 



Hour 


PI 


F2 


P3 


P4 


P5 


P6 


P7 


PS 


P9 


pi a 


Pll 


PI 2 


PI 3 


P14 


PI 5 


PI 6 


PI 7 


SE. 


Load 


1 


173 


173 


173 


173 


173 


130 


130 


130 


130 


105 


105 


105 

















505 


1700 


2 


153 


153 


153 


153 


153 


114 


105 


105 


105 


178 


178 


178 

















475 


1730 


3 


159 


159 


159 


159 


159 


118 











184 


184 


184 


55 


55 


55 


30 


30 


425 


1690 


4 


139 


139 


139 


139 


139 


104 











161 


161 


161 


98.8 


98.8 


98.8 


63 


63 


640 


1700 


5 


15* 


15* 


15* 


156 


156 


116 


105 


105 


105 


105 


105 


105 


55 


55 


55 


30 


30 


365 


1750 


6 


175 


175 


175 


:": 


175 


131 


131 


131 


131 











110 


110 


110 


60 


63 


490 


:s:o 


7 


183 


183 


183 


183 


183 


137 


137 


137 


137 


105 


105 


105 


55 


55 


55 


30 


30 


430 


2000 


S 


173 


173 


173 


173 


173 


129 


129 


129 


129 


201 


201 


201 


110 


110 


110 


60 


60 


540 


2430 


9 


183 


183 


183 


183 


183 


137 


137 


137 


137 


210 


210 


210 


110 


110 


110 


60 


63 


430 


2540 


10 


190 


190 


190 


190 


190 


142 


142 


142 


142 


210 


210 


210 


110 


110 


110 


63 


63 


370 


1603 


11 


199 


199 


199 


'.99 


199 


149 


149 


149 


149 


210 


210 


210 


110 


110 


110 


63 


60 


300 


2670 


12 


189 


189 


189 


189 


189 


141 


141 


141 


141 


210 


210 


210 


110 


110 


110 


60 


60 


380 


2590 


13 


189 


189 


189 


189 


189 


141 


141 


141 


141 


210 


210 


210 


110 


110 


110 


63 


60 


380 


2590 


14 


184 


184 


184 


184 


184 


137 


137 


137 


137 


210 


210 


210 


110 


110 


110 


60 


60 


420 


2550 


:; 


193 


193 


193 


193 


193 


144 


144 


144 


144 


210 


210 


210 


110 


110 


110 


60 


60 


350 


2620 


16 


197 


197 


197 


197 


197 


147 


147 


147 


147 


210 


210 


210 


110 


110 


110 


63 


60 


320 


2653 


17 


203 


203 


203 


203 


203 


152 


152 


152 


152 


210 


210 


210 


110 


110 


110 


63 


63 


270 


2700 


IB 


182 


182 


182 


182 


182 


13* 


13* 


136 


136 


210 


210 


210 


110 


110 


110 


6Q 


60 


440 


2530 


19 


179 


179 


179 


179 


179 


133 


133 


133 


133 


208 


208 


208 


110 


110 


110 


60 


60 


470 


2500 


20 


184 


184 


184 


184 


184 


137 


137 


137 


137 


210 


210 


210 


110 


110 


110 


60 


60 


420 


2550 


21 


190 


190 


190 


190 


190 


142 


142 


142 


142 


210 


210 


210 


110 


110 


110 


60 


60 


370 


7633 


22 


197 


197 


197 


:?" 


197 


147 


147 


147 


147 


210 


210 


210 


110 


110 


110 


60 


60 


320 


2650 


23 


172 


172 


172 


172 


172 


128 


128 


128 


::s 


200 


200 


200 


55 


55 


55 


30 


30 


545 


2200 


24 


191 


191 


191 


191 


191 


143 


143 


143 


143 


105 


13: 


::■: 

















365 


". S43 



It is obvious when the EENS is restricted from 50 MWh to 

spec 

10 MWh the total operating cost is increased from 14870493 
FNR/day to 14946074 INR/day. (i.e). SR al is more when 
EENS is 10 MWh than 50 MWh. 

spec 

VII. Conclusion 

This paper proposed a composite power system network 
equivalent technique to get the solution for reliability con- 
strained unit commitment schedule. The reliability indices 
such as LOLP and EENS are calculated using probabilistic 



approach for the composite power system model incorporat- 
ing load forecast uncertainty. By incorporating transmission 
line outages and load forecast uncertainty an accurate esti- 
mate of spinning reserve requirement is obtained. The re- 
sults found, when compared with those obtained by other 
GA, validate the application of enhanced inherited cross- 
over GA to the reliability constrained unit commitment prob- 
lem. An improvement is observed in the convergence and in 
the quality of the best solution found. The proposed model 
is validated on LEEE RTS 24 bus and a South Indian Tamilnadu 
86 bus system, thereby efficient commitment schedules are 
obtained for the forecasted load. 
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